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^Abstract 
< 

Sweeping is a powerful and versatile method of designing objects. Boundary of volumes (henceforth envelope) obtained by sweeping 
\l solids have been extensively investigated in the past, though, obtaining an accurate paramctrization of the envelope remained 
computationally hard. The present work reports our approach to this problem as well as the important problem of identifying 
HH self-intersections within the envelope. Paramctrization of the envelope is, of course, necessary for its use in most current CAD 
systems. We take the more interesting case when the solid is composed of several faces meeting smoothly. We show that the face 
>^ structure of the envelope mimics locally that of the solid. We adopt the procedural approach at defining the geometry in this work 
C/j which has the advantage of being accurate as well as computationally efficient. The problem of detecting local self- intersections 
O is central to a robust implementation of the solid sweep. This has been addressed by computing a subtle mathematical invariant 
which detects self-intersections, and which is computationally benign and requires only point queries. 
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1. Introduction 

In this paper we focus on the problem of computing an 

accurate paramctrization of the boundary of the volume ob- 
tained by sweeping a solid in along a trajectory and that 
of detecting local self-intersections in the envelope. Sweep- 
ing is an operation of fundamental importance in geomet- 
ric design. It has applications like numerically controlled 
machining verification [3,4,17] and robot motion planning 
[23,24]. There have been several approaches to computing 
the boundary of swept volumes in the past. The works [6,21] 
formulate the problem using the rank deficiency of the Ja- 
cobian, [3 5] compute the envelope by solving sweep differ- 
ential equations, [20] uses inverse-trajectories for deriving 
a point membership test for a point to belong to the enve- 
lope. In [22] the authors give a close approximation of the 
envelope by restricting the trajectories to piecewise screw 
motions. 

Despite the extensive research done in the past in this 
area, computing an accurate paramctrization of the en- 
velope has remained an unsolved problem due to known 
mathematical and computational difficulties [15]. In this 
work we attempt to arrive at an accurate parametrization 
of the envelope through the procedural approach, which is 
an abstract way of defining surfaces and curves when closed 
form formulae are not available. The procedural paradigm 
exploits the fact that from the users point of view, a para- 



metric surface is just a map from to M.'^ and hence can 
be represented in a computer by a procedure which takes 
as input G and returns z) G M.^. Higher or- 

der derivatives of the surface can be returned similarly. The 
definition of splines through the De Casteljau's algorithm 
is an example of procedural parametrization. The authors 
in [14] compute the intersection curve of two parametric 
surfaces by procedural approach. 

The second problem that we tackle in this paper is that 
of detecting local self-intersections. Self-intersections cause 
anomalies in the envelope. There have been rather few at- 
tempts at solving this problem in the past. The paper [2] 
proposes an efficient and robust method of detecting global 
and local self-intersections by checking whether the inverse- 
trajcictory of a point intersects the solid. In the paper [7] 
global and local self-intersections are detected by comput- 
ing intersection of curves of contact at discrete time steps. 
This has the disadvantage of being computationally expen- 
sive. In the paper [18] self-intersections are accurately quan- 
tified and detected but their method is limited to sweep- 
ing tools for NC machining verification. The method em- 
ployed by [25] for detecting local self-intersections is based 
on point set data and could be computationally expensive. 
In [9] the author solves the problem of detecting local self- 
intersections for sweeping planar profiles. In this paper we 
propose a novel test for detecting local self-intersections 
which is based on a subtle mathematical invariant of the 



envelope. It has the advantage of being computationally 
efficient and requires only point queries. 

The paper is organized as follows. In Section 2 we de- 
scribe the input to the sweeping algorithm, in Section 3 
we discuss the overall framework for the computation of 
the envelope, in Sections 4 and 5 we study the mathe- 
matical structure of the envelope and quantify local self- 
intersections, giving a test for detecting them, in Section 6 
we analyse the case when the envelope is free from self- 
intersections. In Section 7 we describe the algorithm for 
computing the procedural parametrization of the envelope. 
We conclude the paper in Section 8. 

2. Preliminaries 

This section outlines the basic representational struc- 
tures associated with the problem. Subsection 2.1 describes 
the boundary representation of a solid which is typical to 
many CAD systems and subsequent sections, the basic in- 
puts and outputs of the sweep algorithm. In Subsection 2.2 
we define the trajectory in along which the solid is swept. 
Next, in subsection 2.3 we define the various mathematical 
sub-entities which make up the envelope. 

2.1. Boundary representation of a solid 

Boundary representation, also known as Brep, is a popu- 
lar and standard method of representing a 'closed' solid M 
by its boundary dM. The boundary dM separates the inte- 
rior of M from the exterior of M. dM is represented using 
a set of faces, edges and vertices. See figure 1 for a Brep of 
a solid where different faces are coloured differently. Faces 
meet in edges and edges meet in vertices. The Brep of a 
solid consists of two interconnected pieces of information, 
viz. the geometric and the topological. 
Geometric information: This consists of geometric enti- 
ties, namely, vertices, edges and faces. A vertex is simply a 
point in M"^. An edge is obtained by restricting the under- 
lying parametric curve by a pair of vertices. A parametric 
curve in is a continuous map 7 : E — >■ E'^ . The curve 7 is 
called regular at sq G E if 7 is differentiable and ^\so 7^ 5- 
Here s is the parameter of the curve. An edge is derived 
from the underlying curve by suitably restricting the pa- 
rameter s to an interval [a, b]. Further, it is required that the 
edge (more precisely, the underlying curve) is regular at all 
points in the interval [a, b] and devoid of self- intersections. 

Similarly, a face is obtained by restricting the underly- 
ing parametric surface by a set of edges. A parametric sur- 
face is a continuous map S : E^ — >■ E'^. The surface S is 
said to be regular at {uq, vq) G E^ if S* is differentiable and 
fil(«o,i'o) G and If |(«o,i>o) G E^ are linearly indepen- 
dent. Here u and v are the parameters of S. A face is de- 
rived from a surface by suitably restricting the parameters 
u and V inside a 'domain'. As expected, it is required that 
the face is regular at all points in the domain and devoid 
of self-intersection. 
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Fig. 1. A Brep of a solid 



Topological information: The topological/combinatorial 
information consists of spatial relationships between dif- 
ferent geometric entities, i.e., the adjacency between faces, 
the incidence relationships between faces and edges and 
so on. In figure 1, for example, the orange and the green 
face are adjacent. Another important component is the 
orientation for each face, that is, a consistent choice of 
outward-normal for that face. The orientation of a regular 
face is a choice of a unit normal from amongst yl"^!"!! and 

— ||f"xg"|| where S is the underlying parametric surface. 
All the faces bounding the solid are oriented so that the 
unit normal at each point on each face is pointing towards 
the exterior of the solid. 

Conceptually, a Brep through its 'global' topological in- 
formation glues the 'local' geometric entities which come 
equipped with associated mathematical parametrizations. 
Note that, the regularity assumptions on the geometric en- 
tities guarantee that the tangent space at every point on an 
edge or a face is of the right dimension. Typically, one also 
imposes higher-order 'parametric' continuity requirements 
which are denoted by C'^ where k refers to the order of con- 
tinuity. For the sake of simplicity, throughout this paper, 
we will assume that the edges and faces bounding the solid 
are regular of class C*^ for some fc > 2, i.e. the underlying 
parametrizations are twice differentiable with continuous 
second order derivatives. Note that, however, these do not 
rule out, e.g., adjacent faces meeting along sharp edges. 

2.2. A trajectory in E^ 

A trajectory in E'^ is a 1-parameter family of rigid mo- 
tions in defined as follows. 

Definition 2.1 A trajectory in M.^ is specified by a map 
h : [0,1] ^ {SO{3),R^),h{t) = {A{t),b{t)) where A{t) G 
5*0(3) 1 , b{t) € E^, A(0) = I, 6(0) = 0. The parameter t m 
this definition represents time. 

For technical convenience, we assume that h is of class C'', 
for some k > 2. 



1 50(3) = {X is a 3 x3 real matrix|X'X = I,det{X) = 1} is the 
special orthogonal group, i.e. the group of rotational transforms. 



2.3. Boundary of the swept volume 

We begin by giving an intuitive description of the bound- 
ary of the swept volume. We will formalize these notions 
in Section 4. Let M be a solid being swept along a given 
trajectory h. By abuse of notation, a point in M will mean 
a point in the interior of M or on the boundary dM of M. 
We denote by Mt the position of M at time t E [0, 1], i.e. 
Mt = {A{t)x + b{t)\x e M}, and by dMt, the boundary 
of Mf. Then Mt is the volume swept by M during 

te[o,i] 

this operation. Our goal is to compute the boundary of this 
swept volume as a Brep, which we will refer to as the en- 
velope. For a fixed point x G Af , consider the trajectory of 
X as the map y : [0, 1] given by y{t) = A{t)x + b{t). 

The trajectory of x describes the motion x in under the 
given trajectory h. Clearly, if x is in the interior of Af , no 
point in the image of the trajectory of x can be on the en- 
velope. Further, at a particular time instant to, only a sub- 
set of points on dMt„ will lie on the envelope. The union of 
such points for all to G [0, 1] gives the final envelope. 

It is clear that, at a given time instant to, only a part of 
(?Af(,-| is in 'contact' with the envelope. To make this more 
clear, fix a point x G dM and the trajectory y of x. The 
derivative of the trajectory of x at a given time instant to, 
that is, ^|to gives the velocity of a; at to- It is easy to show 
that (cf. Section 4) x (more precisely, y{to)) is in contact 
with the envelope at time to only if the velocity of x at to is in 
the 'tangent-space' ofdMtg at y(to). In a generic situation, 
the set of points of dMt^ which are in contact with the 
envelope will be the curve-of-contact. The union of these 
curves-of-contact is called the contact-set or the running 
envelope. Clearly, the total envelope and the contact-set 
are closely related. If all goes well, the envelope is obtained 
from the contact-set by 'capping' it by appropriate parts 
of Mo and Mi, the object at times t = and t = 1. But all 
may not go well. The detection of anomalies is central to 
the use of the algorithm in industrial situations and is an 
important objective of this paper. 

3. Our Approach/Framework 

In this section we briefly describe the overall framework 
for computing the sweep surface i.e. the envelope as a Brep. 
We continue using the notation from the previous section 
where M denotes the Brep/solid and h denotes the trajec- 
tory along which M is swept. 

The naive approach to computing the envelope would be 
to discretize time, i.e, to construct a sequence T = {0 = 
ti, . . . ,tk = 1} and constructing the approximate envelope 
as E' = UjA^t . , the union of the translates. The next step 
would be to construct a smooth version E" of E' above, by 
some fitting operation. However, this approach has several 
issues-(i) computation of E' leads to unstable booleans of 
two very close-by objects, leading to sliver-faces, and (ii) 
the fit of E" to the actual E depends on a dense enough 



choice of T which compounds problem (i) above. There are 

other options, but problems remain. 

In this work, we propose a novel approach based on the 
procedural paradigm (cf [11,14]) which has gained ascen- 
dance in many numerical kernels, e.g., ACIS (cf [12]). 

We now describe the basic architecture for our algorithm. 
For this, we use a running example referred to in figure 1 
and figure 2. The object to be swept is M as in figure 1, 
and the output contact-set is C as shown in figure 2. The 
trajectory is roughly helical with a compounded rotation. 

(i) A natural correspondence between the entities of M 
and the entities ofC. 

Every point p of the envelope comes from a curve 
of contact on Mi, for some t, and therefore belongs 
to some entity of M, i.e., a vertex, edge or face. This 
sets up the correspondence between entities of M and 
those of C. The procedural approach attaches a com- 
mon evaluation method to each such entity. Fig2 il- 
lustrates this correspondence. Faces of C which are 
generated by a particular face of M are shown in same 
colour. Curve-of-contact at time f = is shown im- 
printed on the solid in red. 

Along with the geometric definition of each entity, 
we must also construct the topological data to go with 
it. This data is constructed by observing that there 
is a local homeomorphism between a point on C and 
a suitable point on M. 

(ii) Accurate parametrizations of the geometric entities of 
C with 'time ' as one of the central parameters. 

This is achieved through the procedural paradigm 
in which all key attributes/features of the geometric 
entities are made available through a set of associated 
procedures (cf [11,14]). In our case, these procedures 
are based on Newton-Raphson solvers. This is the 
focus of Section 7. 

(iii) Topological and regularity analysis ofC. 

It is quite common to have a sweeping operation 
in which the resulting envelope/contact-set C self- 
intersects. These self-intersections can be broadly 
classified into global and local self-intersections (see 
figure 6). Once an accurate 'local' parametrization 
(as in step 2) of the faces of C is obtained, in prin- 
ciple, global self-intersections can be detected and 
dealt with by well-known (cf [14]) surface-surface 
intersection solvers. A more subtle mathematical 
issue is that of detecting singularities and local self- 
intersections. This is addressed in Sections 4 and 5. 
We are now in a position to define the scope of this pa- 
per. In this work we describe in detail tasks (ii) and (iii) 
described above, namely, detecting local self-intersections 
and obtaining a procedural parameterization of faces. The 
focus of Section 6 is task (i) in the interesting case when 
M is composed of faces meeting smoothly. Other architec- 
tural aspects and solid-modelling implementation will be 
addressed in a later work. 
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Fig. 2. A solid is swept along a helical trajectory 



4. Mathematical structure of the contact-set 

In this section we will study in detail the mathemati- 
cal structure of the boundary of the volume obtained by 
sweeping the solid M along the trajectory. For simplicity, 
we work with a single parametric surface patch S and anal- 
yse the sweep of S under the trajectory h. As explained 
before, we assume that both S and h are regular of class 
for k > 2, and are devoid of self-intersections. In sec- 
tion 6, we lift the results of this section to the interesting 
case when M is composed of several faces/surfaces meeting 
smoothly. For later use, we introduce the following nota- 
tion: the tanget space to a manifold X at a point p € X 
will be denoted by Txip)- 

We begin with the formal definition of the sweep map. 
Definition 4.1 Given S and h, the sweep is defined as a 
map cr : X [0, 1] — > given by a{u, v, t) = A{t)S{u, v) + 
bit). 

Here u, v are the parameters of S. The position of a 
point S{uo,vo) on surface S at time to will be given by 
a{uQ, vq, to) = A{to)S{uo, vo)+b{to) and the velocity of the 
point S{uo,vo) at time to will be given by ^(^0,^0,^0) = 
If l(«o,«o,to) = A'{to)S{uo,vo) + b'{to), where ' denotes 
derivative with respect to t. If N{uo,Vo) is the unit (out- 
ward) normal to S at {uo,vo), then the unit normal to 
Sto at {uo,vo) is given by TV = A{to) N (uo , vo) where, 
= {A{to)S{u,v) + b(to)\{u,v) e M^} is the position of 
the surface at time instant to- In order to formally define 
the contact-set, we look at the extended sweep in M'* in 
which the fourth dimension is time [2]. 
Definition 4.2 Given S and h, the extended sweep is de- 
fined as a map it : x [0, 1] — > given by d{u, v, t) = 



{a[u,v,t),t). 

Thus, the sweep a is clearly the extended sweep a com- 
posed with the projection map along the t-dimension. De- 
noting partial derivatives using a subscript, we note that 
Uu,cry and CTt are linearly independent for all {u,v,t) G 
X [0, 1] and & is injective. Hence the image of ct is a 3- 
diniensional manifold. We now define the contact-set. 
Definition 4.3 Given S and h, the contact-set is the set 
of points a{uo,vo,to) such that the line {{xo,yo, zo,t) G 
M.'^\(T{uo,vo,to) = {xo,yo, ZQ),t G [0,1]} is tangent to (the 
imageof)aata{uojVo,to) ~ (xq, j/oi -^o, ^o)- We will denote 
the contact-set byC. 

We will refer to the domain of the map a as the parameter 
space and the co-domain as the object space. Consider now 
the following function / : x [0, 1] — >■ K given by 

fiu,v,t)^(viu,v,t),N{u,v,t)) (1) 

Recalling that V{u, v, t) is the velocity of the point S{u, v) 
at time t and N(u, v, t) is the normal to St at (w, v), we look 
at the zero-set of this function in the parameter space. 
Definition 4.4 The funnel J- is defined as the zero-set of 
the function f specified in Eq. 1, i.e.. J- = {(u, v, t) G x 
[0,l]|/(w,z;,t)=0}. 

In other words, if a point p = {uo, vo, to) G T ^ then the 
velocity at the point (t(p) lies in the tangent space of S'tj, 
at g(j)). The following lemma shows that the contact-set is 
precisely the image of the funnel through the sweep map. 
Lemma 4.1 o(fF) = C 

Proof. Skipped here. □ 
Hence S{u, v) 'contributes' a point (namely, A(t)S{u, v)-\- 
b{t)) to the contact-set at time t iff the triplet {u,v,t) 
satisfies the following condition. 
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Fig. 3. Funnel and contact-set 

f{u,v,t) = (v{u,v,t),N{u,v,t) 







(2) 



For a fixed t, Eq. 2 is a system of one equation in two 
variables u and v, hence, in a generic situation, the solution 
will be a curve. 

Eq. 2 can also be looked upon as the rank deficiency 
condition [6] of the Jacobian of the map a defined in 4.1. 
To make this precise, let 



where 



Ja — Uu CTy at (3) 

\iu,v.t) = ^Wff (",w)andcr„|(„,„,t) =A{t)^{u,v) 
and (Jt\(^u,v,t) = V{u,v,t). Observe that regularity of S 
ensures that has rank at least 2. Further, it is easy to 
show that f{u,v,t) is a non-zero scalar multiple of the 
determinant of J^.. Therefore, Eq. 2 is precisely the rank 
deficiency condition of the Jacobian of a. 

Note that, for a point {uo,vo,to) G K^, the Jacobian 
Ja-\{uo,vo,to) is a map from the tangent space to the ambient 
parameter-space at {uq, wq, ^o) to the tangent space to the 
ambient object space at cr(uo, wq, io)- As already noted, if 
(tto, VQ,to) G then Ja\{uo,vo,to) is rank-deficient and maps 
the 3-dimensional ambient tangent space at {uQ,vo,to), 
(surjectively) onto, a 2-dimensional subspace of the ambi- 
ent tangent space at cr(uo, wo,io)- 

The subset of T for a fixed value of t will in general be 
a curve and will be referred to as the pcurve- of- contact at 
time t and its image through a will be a subset of C which 
will be referred to as the curve- of- contact since it is essen- 
tially the set of points on the surface S where S makes tan- 
gential contact with C at time t. The union of such curves- 
of-contact for all t gives the contact-set C. The curve-of- 
contact at t will be denoted by Ct and the pcurve-of-contact 
will be denoted by ct- Fig. 3 schematically illustrates the 
funnel and the contact-set. 

Before proceeding further, we make the following the 
non-degeneracy assumption ^ that: 

Vpe-F,V/|p^ (0,0,0) (4) 

Further, for ease of discussion, we assume that (i) J- is con- 
nected, and (ii) V{u,v,t) e J^,ifu,fv) (0,0) Our analy- 



^ Examples where this does not hold are (i) a cylinder being swept 
along its axis (ii) a planar face being swept in a direction orthogonal 
to its normal. Such cases can be separately and easily handled. 



sis can be easily extended to situations where these further 
assumptions do not hold. 

An important consequence of the assumption 4 is that 

is a 2-dimensional manifold and, hence, Tj^{p) is 2- 
dimensional at all points p G T . Observe that T is also ori- 
entable as V/ provides a continuous non- vanishing normal. 

Thus, T is topologically nice and regular, However, 
quite often, C has 'anomalies' which arise due to self- 
intersections. One of the main contributions of this paper 
is a subtle, efficiently computable mathematical function 
on J- which allows to identify points on F which give rise 
to anomalies in C. 

The key to our analysis ahead, is the restriction of the 
sweep map tr to F . We will abuse the notation, denote this 
restriction, again by a. So, a : T ^ C. Now, fix a point 
p = {u, v,t) G T and let q — a{p) G C. Since det{J„{p)) = 
0, {cr„(p), cry(p), crt(p)} are linearly dependent. As S is reg- 
ular, the set {(Tu{p),c^v{j>)} forms a basis for the tangent 
space to St- Therefore, we must have at ~ la^ +ma^ where 
I and m are well-defined (unique) and are themselves con- 
tinuous functions of u, v and t. 

Clearly, {au{p),ay{p)) is a natural (ordered) 2-frame in 
the object space at point q (recall that, q = a(p)). Fur- 
ther, let X(p) be any ordered continuous 2-frame (basis) 
of the tangent space Tr{p). Note that, this 2-frame is in 
the parameter space and is associated to the point p. Now, 
through CT, more precisely, Jaip), the frame X{p) can be 
transported to another natural 2-framc a(X{p)) in the ob- 
ject space at the point q. The determinant of the linear 
transformation connecting these two natural frames at q, 
namely (i) (ct„(p), ct^(p)) and (ii) a{X{p)) is the key to the 
subsequent analysis. As we show later, this determinant is a 
positive scalar multiple of the continuous function 6 : J- ^ 
M defined as follows. 

eip) = Ifu + mfy - ft (5) 

Here p = {u,v,t) and and ft denote partial deriva- 

tives of the function / w.r.t. u,v and t respectively at p, 
and I and m are as defined before. Note that is easily and 
robustly computed. 

We state an important result which we will prove in the 
coming sections: 

Theorem 4.1 The function 9 is such that (i) 6{p) < 
indicates thatp is a point of local self intersection as defined 
by most authors, (see [2,7]) and (ii) 9(p) = Q is where the 
rank of J„{Tj^{p)) < 2, and finally (Hi) excision of the region 
{p\9(p) < 0} from the funnel J- simplifies the construction 
of the envelope. 

4.1. A particular frame for Tj^ 

Let p = {u,v,t) G J-. In this section, we compute a 
natural 2-frame X{p) in Tj^{p)- Note that, T being the zero 
level-set of the function / defined in Eq. 1, W f\pl.Tj^{p). 
We set /3 = (-/i,,/„,0) ^ and note that /3-LV/. It is 
easy to see that f3 is tangent to the pcurve-of-contact Cj. 
Let a = V/ X /3 = (-/„/«, + f^). Here x is the 
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Fig. 6. Global and local self-intersection 



Fig. 4. A sphere rotating about an axis which is tangent to itself. 
There is typc-1 L.S.I, but no type-2 L.S.I. 

cross-product in M'^. Clearly, the set {a, f3} forms a basis of 

TAp)- 

Figure 3 illustrates the basis {a, /?} schematically. Ob- 
serve that /3 is tangent to the pcurve-of-contact at time t 
and a points towards the 'next' pcurve-of-contact. 

4.2. The determinant connecting the two frames 

We continue with the notation developed earlier. We have 
a = i-ftfu, -ftfv, fu + fv) and fi = (-/„, /„, 0). Hence, 

JcyOl = -ftfu<Ju - ftfyC^v + {fu + fv)<^t 

= (-/Ju + Kfl + fl))<Tu + {-ftf. + m{fl + fl))a, 

Jal3 = —fv<yu + fu'^v 

So, {Jo-a, Jo-/?} can be expressed in terms of {ctujCu} as 
follows 



-ftfu+l{f'a+fv) ~fv 

-ftfv+m{fl + fl) fu 



V{p) 



Note that, 



det{V{p)) = {fl + fl,){lfu + m/„ - /,) (6) 
^{fl + fl)0{p) (7) 

4.3. Singularities of C 

In this subsection, we propose an efhcient test for detect- 
ing singularities on the contact-set. See Fig. 4 for an exam- 
ple. Clearly, the detection of singularities is important in 
practice. 

We start with the following definition. 
Definition 4.5 We say that the sweep causes a singularity 
if the composite map T^C^R^ fails to be an immersion. 
In other words, the sweep causes a singularity if 3p e J" 
such that the rank of Ja{p){Tj^{p)) is less than 2. Following 
the standard usage [8] , in this case we say that the point p 
is a critical point. 

Lemma 4.2 A point r Cz J- is a critical point iff 0{r) = 
iff the rank of J„{r){Tj^{r)) is less than 2. 



Proof. By equation 7, we have det(I?(r)) = iff 9{r) = 0. 
Recall that, as shown earlier, {a(r),/3(r)} is a basis of 
Tj^{r), and {au{r),ay{r)} is also a basis of the 2-frame 
associated at a{r). As I'(r) is the matrix expressing 
X = {J^{r){a{r)),J„{r){a{r)) in terms of {(T„(r), (J„(r)}, 
det(I?(r)) = iff rank of J<.(r)(7>(r)) is less than 2. Thus, 
r is a critical point iff 9{r) = iff the rank of Ja{i'){Tr{r)) 
is less than 2. □ 
Note that the above lemma proves part (ii) of theo- 
rem 4.1. 

Lemma 4.3 A sweep causes a singularity if there exists 
points p and q on F such that 9(jp) < and 6{q) > 0. 
Proof. As is a continuous function on T, the existence of 
p and qonT with the required properties implies existence 
of another point r ^ F such that 9{r) ~ 0. By the previous 
lemma, this implies that the sweep causes a singularity. □ 

The above lemma leads to a computationally efficient test 
for detecting singularities: namely, evaluating 9 at sampled 
points on T and checking if it changes sign on !F. 

The analysis done so far helps us detect singularities on 
the contact-set. In the next section we will perform a de- 
tailed analysis of local self-intersections which is topologi- 
cal in nature. Towards this, note that all points in the non- 
critical set may not lead to points on the envelope of the 
swept volume. For some p = (uo,i'o,io) G -^i cr(p) may lie 
in the interior of the solid Mt of which the surface patch St 
is a part of, for some t in neighbourhood of t^ . In that case 
<t{jp) will not be on the envelope. Fig. 5 shows two sweep- 
ing examples with self-intersections. Curves-of-contact at a 
few time instances are shown. In the next section we focus 
on identifying such points. 

5. Topological and regularity analysis of C 

Quite often, the anomalies on C arise due to self- 
intersections. If C has self- intersections, it needs to be 
trimmed to obtain the envelope of the swept volume [2]. 
Self-intersections can be broadly be classified into global 
and local. Fig. 6 illustrates the difference between global 
and local self-intersections schematically. 

If there are only global self- intersections occurring on C, 
the composite map A C ^ fails to be an injection. 
However, it is an immersion (see [8]), i.e. Vp € JF , the rank 
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Fig. 5. Examples of local self-intersection showing curves-of-contact at few time instances: (a) A cylinder with blended edges undergoing 
translation and rotation about i-axis (b) An ellipsoid undergoing translation along a curvilinear path 



of Ja{p){Tr{p)) is 2. In principle, global self-intersections 
can be detected by surface-surface intersection, (see [14]) 

The case of local self-intersection is more subtle as it 
leads to singularities on C. Clearly, the detection of local 
self-intersections is also central to a robust implementation 
of the solid sweep in CAD systems. 

In literature [2,7], local self-intersections have been quan- 
tified by looking at points in the contact-set which lie in the 
interior of the solid Mt for some time instant t. Clearly, such 
a point cannot be on the envelope of the swept volume. This 
approach was used in [2] for detecting local and global self- 
intersections, where the authors used implicit representa- 
tion of the surface bounding the solid which is being swept. 
We adapt this concept to when the surface of the solid is 
represented parametrically. We will refer to this type of lo- 
cal self-intersection as type-2 L.S.I.. It turns out that type- 
2 L.S.I, is intimately related to the analysis carried out in 
the previous section. To make this connection precise, we 
first introduce another type of local self-intersection. For 
lack of a better name, this is called as type-1 L.S.I. 

5.1. Type-1 local self-intersections 

Definition 5.1 A type-1 L.S.I, is said to occur at a point 
p&T ifd{p) < 0. 

Thus, Type-1 L.S.I is our classification of a local self 
intersection. We will see in subsection 5.3 that a for a Type- 
1 L.S.I point p, the image a{p) does not lie on the envelope 
of the swept volume. 

5.2. Type-2 local self-intersection 

In order to define type-2 L.S.I, we first describe the in- 
verse trajectory corresponding to a given trajectory [2,20]. 

Given a trajectory as in definition 2.1 and a fixed point x 
in object-space, we would like to compute the set of points 
in the object-space which get mapped to x at some time 
instant. This set can be computed through the inverse tra- 
jectory defined as follows. 



Definition 5.2 Given a trajectory h, the inverse trajec- 
tory h is defined as the map h : [0, 1] — > (50(3), M'^) given 
byh{t) = {A\t),-A\t)b{t)). 

Thus, for a fixed point x e M'^, the inverse trajectory of 
X is the map y : [0,1] — )■ given by y{t) = A*{t){x — 
b{t)). The range of y is {A\t)x - A^{t)b{t)\t e [0, 1]} = 
{z e R^\3t e [0, l\,A{t)z + b{t) = x}. We will denote the 
trajectory of a; by y : [0, 1] y{t) = A{t)x + b{t). We 

now note a few useful facts about inverse trajectory of x. 
We assume without loss of generality that ^(^o) — I and 
b{to) = 0. Denoting the derivative with respect to t by ', we 
have 

m^A\t){x-b{t))-A\t)b{t) (8) 
Since A G SO{3) we have, 

A\t)A{t) = lyt (9) 
Differentiating Eq. 9 w.r.t. t we get 

A\t)A{t) + A\t)A{t) = oyt (10) 
i*(io) + A(to) = (11) 
Differentiating Eq. 10 w.r.t. t we get 

A\t)Ait) + 2A\t)A{t) + A\t)A{t) = oyt 

A\to)-^2A\to)Aito)+A{to)^0 (12) 

Using Eq. 8 and Eq. 11 we get 

§{ta) = -A{to)x - b{to) - -y(io) (13) 

Differentiating Eq. 8 w.r.t. time we get 

m = A\t){x - bit)) - 2A\t)b{t) - A\t)b{t) (14) 

Using Equations 14, 11 and 12 we get 

y{to) ^ -y{to) + 2A{to)y{t„) (15) 

We now define type-2 L.S.I. The surface S is the boundary 
of the solid M being swept. We will refer to the interior of 
M by Int{M) and exterior of M by Ext{M). 
Definition 5.3 A type-2 L.S.I is said to occur at a point 
(uq, Vo,ta) if the inverse trajectory of the point cr(uo, "^01^0) 
intersects Int{Mtg) .(see [2]) 
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Fig. 7. Type-2 local self-intersection 

Fig. 7 illustrates type-2 L.S.I, schematically where y is 
the inverse trajectory of the point x G St„ and tt is the 
projectionof yon 5(0. Suppose the point a; = cr(uo,uo,io) G 
C. Let \{t) be the signed distance of y{t) from surface Sto- 
If the point y{t) is in Int{Mtg), Ext{Mtg) or on the surface 
Stg , then A(t) is negative, positive or zero respectively. Then 
wehavey(t) — 7r(i) = X{t)N{t), where ^{t) is the projection 
of y(t) on Stg along the unit outward pointing normal N(t) 
to Stg at TT{t). Then, the following relation holds for A. 

X{t) = m-At),Nit)) (16) 

We now give a necessary and sufficient condition for type- 
2 L.S.I, to occur. 

Lemma 5.1 Type-2 L.S.I, occurs at a point p = {uQ,vo,to) G 
J" if and only if either of the following conditions hold 

(i) X{to) = (^-a + 2Av,N'j + kv'^ < where K is the 

normal curvature of Stg at{uo,VQ) along velocity V{p), 
N is the unit length outward pointing normal to Stg 
at {uo, vq) and v'^ = {V{p), V{p)) . 
(ii) X{t) is negative for some t in some nbd of to. 
Remark: The statement of the above lemma (except for 
the insightful expression of X(to)) is similar in spirit to the 
key Theorem-2 in [2] in the context of implicitly defined 
sohds. 

Proof. Differentiating Eq. 16 with respect to time and de- 
noting derivative w.r.t. t by ', we get 

X{t) = {g{t) - ir{t), N{t)) + (y(i) - 7T{t), N{t)) (17) 

x{t) = im ~ m,N{t)) + 2 (^(i) - m,N{t)) 
+ (^m-At),N{t)) (18) 

At t = to, y{to) = 7r(to). Since y{to) - V{p)LN{p), it 
follows from Eq. 13 that y{to)-\-N{p). It is easy to verify 
that 'k{to) — y{to)- Hence, 

A(to) = A(io) - (19) 

From Eq. 18 and Eq. 15 it follows that 

X{to) = {^{to)-n{to),N{to)) 

= (-y(io) + 2A{to)y{to) - n{to),Nito)) (20) 

Since 7r(t) G Stg^t in some neighbourhood U of to, we 
have that {7r{t),N{t)) = 0,Vi £ U. Hence {TT{t),N{t)) + 

^7r(t),iV(t)^ = 0,V< e U. Hence - {f[{to), N(to)) = 



(^n{to),Nito)) = {n{to),g*{n{to))) = {y{to),g* (yito))) = 

{V{p),G* {V{p))) = Kv'^. Here G*{y) is the differential of 
the Gauss map, i.e. the curvature tensor of Stg at point 
X. Using this in Eq. 20 and the fact that y{to) — o-(p), 
y(^o) = <^{p) we get 

A(io) = {-aip) + 2A{to)V{p),N{to)) + KV^ (21) 

From Eq. 19 and Eq. 21 we conclude that if X{to) < the 
point X = y{to) — a{uo,vo,to) is a local maxima of the 
function A and the inverse trajectory of x intersects with 
interior of the solid Alf^ causing type-2 L.S.I. Similarly, if 
X{to) > we conclude that x is a local minima of A and the 
inverse trajectory of x does not intersect with the interior 
of Mtg and there is no L.S.I, occurring at x. 

However, if X{to) is zero, one needs to inspect a small 
neighbourhood of tg to see if 3t such that X{t) < in order 
to check for type-2 L.S.I. □ 

If A = at a point, the structure of the contact-set C is 
unknown at that point. We will see in the next subsection 
that at such a point, C has singularity. 

5.3. Relation between type-1 L.S.I, and type-2 L.S.I. 

In this subsection we will see that type-2 L.S.I, implies 
type-1 L.S.I, at any point p = (mq, vo, to) e J- 
Lemma 5.2 8{p) = X{to). 
Proof. Recalling definition of 9{p) from Eq. 5 

Ifu + mf, - ft^ (iNu + mN,,v) + (n , IVu + mV^) 

-(^Nt,v)-(N,Vt) 

Here Nu = Q*{o-u) and Ny = Q*{av) where Q* is the shape 
operator (differential of the Gauss map) of Sto at (wo,fo)- 
Also, Vu = AtSu and Vy = AtSy. Assume without loss of 
generality that A{to) = / and h{to) = 0. Using Eq.ll and 
the fact that at = IcTu + tict^ we get 

Ifu + mf, -ft = {g*V, V) + 2 {AtV, N) - {Vt, N) 

^ Kv"^ + {2AtV -Vt,N) (22) 

From Eqs. 21 and 22 and the fact that ^ = Vt we get 

e{p) - Ifu + mf,-ft = x a 

Let p be such that 9{p) < 0. By lemma 5.2, A(to) < 0. 
Further, by lemma 5.1 a type-2 L.S.I, occurs at p. This 
proves parts (i) and (iii) , and hence completes the proof of 
theorem 4.1. 

Fig. 8 schematically illustrates the region on the funnel 
where local self-intersection occurs and the corresponding 
region on the contact-set. A curve-of-contact and the cor- 
responding pcurve-of-contact is shown in red colour at a 
time instant to when a local self-intersection occurs. The 
shaded region there corresponds to 9 < 0. Thus sign{9) 
changes (from — ve to -l-ve) as one moves from the interior 
to the exterior of the shaded region. Of course 6* = A = 
on the boundary of the shaded region where C has a singu- 
larity. 
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Fig. 8. Region of local sclf-intcrsection 



Fig. 4 shows an example of sweeping which clearly 
demonstrates the subtle difference between type-1 L.S.I, 
and type-2 L.S.I. In Fig. 4 a sphere is rotated about an 
axis which is tangent to itself. The contact-set is shown 
in green and the curve-of-contacts on the contact-set at 
initial position {t = 0) and final position (t = 1) are also 
shown in blue. In Fig. 4 there is type-1 L.S.I, but no type-2 
L.S.I, since no point of the contact-set intersects with the 
interior of solid at any time. 

We now present two examples of sweeping which result in 
local self-intersections. We apply the test for type-1 L.S.I, 
on these examples to demonstrate its effectiveness. 
Example 5.1 Consider a solid cylinder along the y-axis 
with height 2.5 and radius 2, parameterized as S{u,v) = 
(2cosw,M, — 2sinw), u g [—1.25,1.25], v £ [—TT,Tr] being 
swept along the trajectory given by h{t) = {A{t),b{t)) 

1 



where A{t) 



3 cos {^t 



cos(0.l7ri) -sin(0.l7ri) 
sin(0.l7rt) cos(0.l7ri) 



and b{t) 



3 sin 







t € [0, 1]. The resulting envelope has lo- 



cal self-intersections as illustrated in Fig. 5(a) which shows 
the cylinder with blended edges and curves-of-contact at 
few time instances. Type-1 L.S.I is detected at time t — 0.1 
by the test given in lemma 5.1 at point p = (u = 0.18, v = 
1.53, t = 0.1). e{p) = -2.378. 

Example 5.2 Consider an ellipsoid with axes lengths 3, 1, 
and \ parameterized as S{u,v) = (— 3cos(u) cos(w), sin(u), 
cos(u) sin(ti)), u G ["fif]?"*^ G [^^: being swept along 
the curvilinear trajectory given by h{t) — {A{t),b(t)) 



where A(t) = Tit and b{t) = 



3 cos 
3 sin 







The curves-of-contact are shown in Fig. 
L.S.I, is detected at time t ~ 0.8. 9{p) 
(-0.791,-0.157,0.8). 



t e [0,1]. 



5(b). Type-1 
= -11.864 at 



6. Mathematical structure of the smooth case 

In this section, we consider the smooth case where the 
solid M is composed of faces meeting smoothly. As usual, 
each face (or the associated surface patch) is smooth (of 
class C'^ for A: > 2). Further, adjacent faces meet smoothly 
at the common edge. This is referred to as continuity [13] 
which formally means that the unit outward normals to the 
adjacent faces match on the common edge. Similarly, at a 
vertex, all the unit outward normals to faces incident on 
this vertex are identical. The solid shown in figure 1 is such 
a solid. 

Consider a sweep of AI along a trajectory h which causes 
no self- intersections/anomalies on the contact-set C. As de- 
scribed in Section 3, every point p on C comes from a curve 
of contact on M and therefore is associated to a point q 
of M. Let TT : C AI [p ^ tt{p) — q) denote this natu- 
ral map. For every p e C, 7r(p) belongs to some geometric 
entity of M, i.e., a vertex, edge or face. This sets up the 
natural correspondence between geometric entities of C and 
that of M . For a face F of M, let Cp denote the part of C 
which corresponds to the face F under this correspondence. 
For example, in figure 2, the green face on the solid M cor- 
responds to multiple green faces of C. Clearly, the map tt 
restricts naturally from Cp ^ F . 

There are situations in which, for example, an edge (or 
a part of it) on M remains on the boundary for a while 
and thus, 'sweeps' a face on C. For simplicity, we assume 
that such situations are ruled out. In other words, no lower 
dimensional geometric entity of AI gives rise to a higher 
dimensional geometric entity on C. This is the case for the 
sweep operation illustrated in figure 2. 

6.1. Local similarity within a face 

Firstly, recall that each face F of M is derived from an 
underlying surface Sp by restricting the parameters of Sp 
to a suitable domain Dp. Now, by applying the 'local' anal- 
ysis of Section 4 to the surface Sp,we have 
Lemma 6.1 The setCp has no self-intersections and is a 
smooth manifold. Further, let p G Cp correspond to q € F 
at time t. Then, the unit normal N{p) to Cp at p is simply 
A{t)N{q) where N{q) is the unit normal to F at q. 

Further, we would like to show that 'locally', Cp has the 
same topology as that of F. More precisely, the natural 
map TT : Cp F is a, local homeomorphism onto its image. 
Thanks to the 'local' nature, we may analyse this via the 
underlying surface Sp. For ease of notation, we sometimes 
omit the reference to F and freely use notations from Sec- 
tion 4. As the sweep a is free of self-intersections, the key 
map a : T Cp is a bijective immersion (recall that, in 
Section 4, F is the funnel defined via the function / = 0). 
Thus, by the inverse function theorem, it is invertible via 
a continuous inverse. Therefore, in order to show that tt : 
Cp F is SL local homeomorphism onto its image it suffices 
to prove the following lemma. 



n 



Lemma 6.2 The natural map tt' : J" — >• F defined as: for 
p = {u,v,t) e 7r'(p) = S{u,v), is a local homeomor- 

phism. 

Remark The map tt' is simply the composition of a and tt. 
Proof. Let p = (wq, '-"oj ^o) G Here, the assumption that 
/t ^ at p is very crucial. Firstly, by the implicit function 
theorem, there exists a neighbourhood O of {uo,vo) and 
a continuous function t = g{u,v) defined on O such that 
y(u,v) e O, f{u,v,g{u,v)) = 0. Further, the set 0{p) = 
v,g(u, v)) I (u, v) G O} is a neighbourhood of p in J^. 
On this neighbomhood of p, the function tt' is a bijcction 
and invertible. This easily follows from the fact that, for 
every (m, v) G O, there is a unique time t, namely, t = 
g{u, v), such that {u, v, g{u, v)) e 0{p). □ 

6.2. Local similarity across faces 

We begin by studying the variation of the unit normal 
across faces of C. Let Ct denote the curve of contact of C at 
time t. Let p E Ct he such that p is common to (only) Cpi 
and where Fi and F2 are two distinct faces of M. 
Lemma 6.3 The faces Fi and F2 are adjacent in M. Fur- 
ther, the normal to Cp^ atp is identical to the normal to 
at p. 

Proof. Suppose p corresponds to q. Clearly, q is common 
to (only) Fi and F2. Thus Fi and F2 are adjacent in M. 
Further, by continuity, the normal to Fi at q is identical 
to the normal to F2 at q. By lemma 6.1, it is clear that, the 
normals to , and to , at p are identical. □ 

Thus, the adjacencies on C are the 'same' as the ad- 
jacencies on M. See figure 2, which effectively illustrates 
this through colours. Further, the adjacent entities of C 
meet smoothly across common lower-dimensional entities. 
In other words, C is also of class . Recall that the overall 
envelope may be obtained from the contact-set C by sim- 
ply capping the appropriate parts of the solid M at the ini- 
tial and final position. Therefore, the topological structure 
of the contact-set and hence, that of the envelope, mimics 
that of the solid. 

The following theorem summarizes the analysis so far. 
Theorem 6.1 The map w from the contact-set/envelope to 
the solid is an adjacency-respecting local homeomorphism 
onto its image. 



and we get a local parameterization of C in A/" by il>{u, v) = 
u(u,v,t(u,v)). Tc{p) is spanned by %Iju = Ou + (^ttu and 
i'v = + cTttv Since p G F \>y lemma 4.1, ot is in the 
space spanned by (t„ and (t„. Let at = lau + may. Hence we 
express basis for Tc{p) in terms of basis of Ts{q) as follows 



where M = 



mtu 1 -l- mty 



M 



(23) 



. The unit normal to C is 



given by N{u,v) = A{t{u,v))N{u,v) = N{u, t!) where N 
is the unit normal to St- Hence, Nu = and Ny = Ny. 
Further, 



On CT'L 



w 



(24) 



where W is the Weingarten matrix of S at point q whose 
determinant gives the Gaussian curvature of S at q{see [1]). 
From Eq. 23 and Eq. 24 we have 



N,, AT,, 



(25) 



So the Weingarten matrix of C with respect to parame- 
terization is given by M~-'^W and the Gaussian curva- 
ture is given by dei(M~^W). From Eq. 2 and the fact that 
Nt = Vu = Vy=Owe note that / = (TV, V), /„ = (7V„, V), 
f, = {Ny,V) a.ndft = {N,Vt).So, 



det{M-^] 



1 



l + ltu + rntv ft - Ifu - mfy 
^ {N,Vt) 

{N,Vt)-{lNu + mNy,V) 

^ {N,Vt) 

(iV,T4)-(W(T/),F) 

The Gaussian curvature of C is computed as 

{N,Vt) 



dei(M-iW) = 



det{W) (26) 



{N,Vt)-{W{V),V) 

where, (iV, Vt) is the curvature of the trajectory scaled by 
IjVfll, (yV{V),V) is the normal curvature of 5 at p in di- 
rection V scaled by HV^H^ and det(W) is the Gaussian cur- 
vature of S at p. 

7. Envelope computation 



6.3. Curvature ofC 

In the special case when the trajectory h consists only of 
translations, i.e. A(t) = I yt, the Gaussian curvature of C 
can be expressed in terms of the Gaussian curvature of S 
and the curvature of h. Since A{t) = iyt,au = Sy and = 
Sy. Also, Vu = Vy = 0. Consider a point p — (uq, Voj^o) € 
F. By definition of F, we have that /(mq, vq, to) = 0. Given 
that V/lp ^ suppose without loss of generality that 
ft\p 7^ 0- Then by the implicit function theorem there ex- 
ists a neighbourhood N oiq= (wq, Vq) such that V(u, v) G 
jV, f{u,v,t{u,v)) = 0. Hence, = and ty = 



In this section we describe the construction of the enve- 
lope £ assuming that it is free from self-intersections and 
hence regular. We obtain a procedural parametrization of 
£. The procedural paradigm is an abstract way of defin- 
ing curves and surfaces. It relies on the fact that from the 
user's point of view, a parametric surface (curve) in is 
a map from M^(]R) to and hence is merely a set of pro- 
grams which allow the user to query the key attributes 
of the surface (curve), e.g. its domain and to evaluate the 
surface(curve) and its derivatives at the given parameter 
value. The procedural approach to defining geometry is es- 
pecially useful when closed form formulae are not available 



for the parametrization map and one must resort to itera- 
tive numerical methods. We use the Newton-Raphson(NR) 
method for this purpose. As an example, the parametriza- 
tion of the intersection curve of two surfaces is computed 
procedurally in [14]. As we will see, this approach has the 
advantage of being computationally efficient as well as ac- 
curate. For a detailed discussion on the procedural frame- 
work, see [11]. 

The computational framework is as follows. For the in- 
put parametric surface S and trajectory h, an approximate 
envelope is first compiitcd, which we will refer to as the seed 
surface. Now, when the user wishes to evaluate the actual 
envelope or its derivative at some parameter value, a NR 
method will be started with seed obtained from the seed 
surface. The NR method will converge, \ipto the required 
tolerance, to the required point on the envelope, or to its 
derivative, as required. Here, the precision of the evalua- 
tion is only restricted by the finite precision of the com- 
puter and hence is accurate. It has the advantage that if a 
tighter degree of tolerance is required while evaluation of 
the surface or its derivative, the seed surface does not need 
to be recomputed. Thus, for the procedural definition of 
the envelope we need the following: 

(i) a NR formulation for computing points on £ and its 
derivatives, which we describe in subsection 7.1 

(ii) Seed surface for seeding the NR procedure, which we 
describe in subsection 7.2 

Recall that by the non-degeneracy assmnption, £ is the 
union of Ct,^t. This suggests a natural parametrization 
of £ in which one of the surface parameters is time t. We 
will call the other parameter p and denote the seed surface 
by 7 which is a map from the parameter space of £ to 
the parameter space of tr, i.e. ^{p,t) = {u{p,t),v{p,t),t) 
and while the point a{'j{p,t)) may not belong to £, it is 
close to £. In other words, j{p,t) is close to J^. We call 
the image of the seed surface through the sweep map a as 
the approximate envelope and denote it by £, i.e. £{p, t) = 
cT{-y{p, t)). We make the following assumption about £. 
Assumption 1 At every point on the iso-t curve of£, the 
normal plane to the iso-t curve intersects the iso-t curve of 
£ in exactly one point. 

Note that this is not a very strong assumption and holds 
true in practice even with rather sparse sampling of points 
for the seed surface. We now describe the Newton-Raphson 
formulation for evaluating points on £ and its derivatives 
at a given parameter value. 



So, given {po, to), we have one equation in two unknowns, 
viz. u(po,^o) and v{pQ,t()). £{po,to) is defined as the in- 
tersection of the plane normal to the iso-t (for t — to) 
curve of £ at £{po,to) with the iso-i(for t — to) curve of 
£ which is nothing but Cj^. Recall that Ctg is given by 
<j{u{p, to), v(p, to), to) where u, v, t obey Eq. 27. Henceforth, 
we will suppress the notation that u,v,u and v arc func- 
tions of p and t. Also, all the evaluations will be understood 
to be done at parameter values {po, to) - The tangent to iso-t 
curve of £ at (^05*0) is given by 



(28) 



d£ da du da dv 
dp du dp dv dp 

Hence, £{po,to) is the solution of simultaneous system of 
equations 27 and 29 



a{u,v,to) - a{u,v,to), 



d£ 
dp 







(29) 



Eq. 27 and Eq. 29 give us a system of two equations in two 

unknowns, u and v and hence can be put into NR framework 
by computing their first order derivatives w.r.t u and v. For 
any given paramc;ter value; (po, to), we seed the NR method 
with the point {u{po,to),v{po,to)) and solve Eq. 27 and 
Eq. 29 for {u{po,to),v{po,to)) and compute £{po,to)- 

Having computed £{p,t) we now compute first order 
derivatives of £ assuming that they exist. In order to com- 
pute we differentiate Eq. 27 and Eq. 29 w.r.t. p to ob- 
tain 



dNdu ^ dNdv 
du dp dv dp' i 



. ^ dV du dV dv 
' du dp dv dp 



' da du da dv 
, du dp dv dp 



da du da dv d£ \ 
du dp dv dp' dp I 

d'^B" 

a{u,v,to) - a{u,v,to), ) = 



= 
(30) 

(31) 



Eq. 30 and Eq. 31 give a system of two equations in two un- 
knowns, viz., 1^ and |^ and can be put into NR framework 

by computing first order derivatives w.r.t. |^ and Note 

that Eq. 30 and Eq. 31 also involve u and v whose compu- 
tation we have already described. After computing |^ and 



can be computed as f" 



dcr dv dS 



&p' dj> v.v^i"i^Ln.c^ 0.0 aj, T a^^g^- St can similarly 

be computed by differentiating Eq. 27 and Eq. 29 w.r.t. t. 

Higher order derivatives can be computed in a similar 
manner. 



7.1. NR formulation for faces of £ 



7.2. Computation of seed surface 



Recall that the points on £ were characterized by the 
tangency condition given in Eq. 2. Introducing the param- 
eters {p, t) of £, we rewrite Eq. 2 V(po, ^o): 

f{u{po,to),v{j)o,to),to) = (N{u{po,to),v{po,to),to), 



V{u{po,to),v{po,to),to)) =0 

(27) 



The seed surface is constructed by sampling a few points 
on the envelope and fitting a tensor product B-spline sur- 
face through these points. For this, we first sample a few 
time instants, say, T = {ti,t2, ■ ■ ■ , t„} from the time inter- 
val of the sweep. For each ti S T, we sample a few points on 
the curve-of-contact Ct . . For this, we begin with one point 
p on and compute the tangent to Cj^ at p, call it Td {p)- 



p + Tci (p) is used as a seed in Newton- Raphson method to 
obtain the next point on Ct^ and this process is repeated. 

While we do not know of any structured way of choosing 
the number of sampled points, in practice even a small 
number of points suffice to ensure that the Assumption 1 
is valid. 

7.3. NR formulation for edges and vertices of£ 

An NR formulation for edges and vertices of £ can be 
obtained in a manner similar to that for faces of £ which we 
described in subsection 7.1. In order to obtain a procedural 
parametrization for edges of £, again seed curves need to 
be computed. 

8. Discussion 

In this work, we have proposed a novel computation- 
ally efficient test for detecting anomalies on the envelope. 
This has been achieved through a delicate mathematical 
analysis of an 'invariant'. We have provided a rich proce- 
dural framework for computing the Brep of the envelope 
along with its accurate parametrization. Another contri- 
bution is a natural correspondence between the gcomet- 
ric/topological entities of the Brep envelope and that of 
the Brep solid. This framework has been implemented us- 
ing the ACIS kernel [12] and has been used to produce the 
running examples of this paper. 

Ongoing work includes, for example, extending the pro- 
posed procedural framework to handle (i) deeper topolog- 
ical information of the Brep envelope, (ii) swept edges (an 
edge on the solid sweeping a face on the envelope) and so on 
(iii) faces meeting with continuity (i.e. sharp edges). We 
also to plan to extend the detection of anomalies to above 
settings and further, trim the appropriate part to obtain 
the final envelope. 

Another exciting future direction would be to analyse 
sweeps in which some numerical invariant associated with 
a a ciirve of contact, varies over time. For example, one 
may imagine sweeping a torus along a trajectory where the 
number of components of the curve of contact changes over 
time. In such a case, one would like to efficiently compute 
deeper topological invariants, say genus, of the envelope. 
Our mathematical analysis coupled with a Morse-theoretic 
analysis appears promising. 
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